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Abstract. This paper describes a conceptual framework for understanding kinetic 
plasma turbulence as a generalized form of energy cascade in phase space. It is 
emphasized that conversion of turbulent energy into thermodynamic heat is only 
achievable in the presence of some (however small) degree of collisionality. The 
smallness of the collision rate is compensated by the emergence of small-scale structure 
in the velocity space. For gyrokinetic turbulence, a nonlinear perpendicular phase 
mixing mechanism is identified and described as a turbulent cascade of entropy 
fluctuations simultaneously occurring at spatial scales smaller than the ion gyroscale 
and in velocity space. Scaling relations for the resulting fluctuation spectra are derived. 
An estimate for the collisional cutoff is provided. The importance of adequately 
modeling and resolving collisions in gyrokinetic simulations is bicfly discussed, as well 
as the relevance of these results to understanding the dissipation-range turbulence in 
the solar wind and the electrostatic microturbulence in fusion plasmas. 
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1. Turbulence: the Symptoms and the Cause 

What is turbulence? Modulo many definitional and interpretational subtleties [T8l I51] . 
turbulence is multiscale disorder, we tend to say that we are dealing with a turbulent 
system if we have detected (measured, observed, simulated, intuited) chaotic fluctuations 
of some field(s) over a broad range of scales. In plasmas, these fluctuating fields 
are the electric and magnetic fields and the distribution function of the particles 
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(either measured directly or accessible partially via its moments: density, flow velocity, 
temperature). So turbulence is defined as a syndrome jl6]: it is identifed by its 
symptoms!) The next logical step is to ask what causes the development of the problem 
in the first place. The short answer is energy injection: in physical systems, turbulence 
is stirred up by some source of energy, which is system-specific and can be in the form of 
direct mechanical forcing (spoon in a tea cup, supernovae in the interstellar medium), 
boundary conditions (airplane wing), or various instabilities feeding on background 
equilibrium gradients (tokamak microturbulence, solar convection, magnetorotational 
turbulence in accretion discs). The fluctuation energy injected into the system is nearly 
always dissipated into heat. Because the dissipation mechanisms available to the system 
have to do with its material properties (microphysics) and are usually unrelated to 
the energy- injection mechanism (macrophysics), there is more often than not a scale 
separation between the energy-injection, or energy-containing, scale (the outer scale) 
and the much smaller dissipation scale (the inner scale). In order to dissipate energy, 
the system has to bridge this gap and one way for this to happen is for the nonlinear 
interactions to fill the intermediate scale range with fluctuations - - giving rise to 
multiscale disorder, or turbulence (there are, of course, other ways, e.g., shock or current- 
sheet formation, but we will not consider them here). 

The simplest illustration of the argument made above is the case of a Navier- Stokes 
neutral fluid, whose velocity field u satisfies 

d t u + u ■ Vu = - Vp + vV 2 u + /, V-u = 0, (1) 

where p is pressure, v the molecular viscosity of the fluid, and the body force / stands 
in for the outer-scale energy injection. The kinetic energy of the fluid then satisfies 

d f d 3 r u 2 f d 3 r ._ l9 , . 

s-v / — |Vw| 2 , (2) 



dt J V 2 J V 

where V is the system volume and s = (1/V) f d 3 r u ■ f is the injected power per unit 
volume. In a stationary state, the injection and dissipation terms on the right-hand 
side of this equation must balance, even though e is finite and viscosity is small, or, 
more precisely, the viscous term in ([1]) is negligible at the outer scale. The balance is 
accomplished by transferring kinetic energy to small scales, where the velocity gradients 
are large, compensating for the viscosity's smallness. The viscous (inner) scale to 
which the energy has to travel in order to be dissipated is, on dimensional grounds, 
L ~ (^/e) 1 / 4 ~ LRe~ 3//4 , where L is the outer scale and Re = u rras Ljv is the Reynolds 
number. The system becomes turbulent when Re 3> 1, i.e., l v -C L, so fluctuations arise 
over a broad band of scalesjfl 

2 We thank T A Yousef for bringing to our attention this analogy, which is particularly apt in fusion 
contexts, where turbulence is indeed a disease that gives rise to anomalous transport, prevents plasma 
confinement and thus hampers humanity's progress toward the hydrogen-powered future. 

3 A one-paragraph review of the Kolmogorov-Obukhov 1941 turbulence theory [32, 38 : If it can be 
assumed (by no means an automatic certainty!) that the energy is transported locally from scale to scale 
[40j , the energy flux through the intermediate scales L > A > i„ (the inertial range) must be constant 
and equal to e. Assuming that fluctuations are isotropic in the inertial range, we have 5u\/t\ ~ e, 
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2. Plasma Turbulence: Entropy, Heating and the Kinetic Cascade 

Can this argument be generalized to plasma turbulence? If the plasma is sufficiently 
collisional, its dynamics is described by a set of fluid equations with diffusive dissipation 
terms [9]. While things become more complicated than for the Navier-Stokes equation 
(multiple fields and species, different diffusion coefficients perpendicular and parallel 
to the magnetic field, interplay between waves and nonlinear interactions), the basic 
principle remains the same: small-scale spatial structure is generated so that the energy 
injected at the outer scale can be transferred to the smaller dissipative scales and 
converted into heat. All this, however, is only valid for fluctuations whose characteristic 
spatial and temporal scales remain collisional, namely fc||A m f p <C I and u <C va, where 
fcii is the typical wavenumber parallel to the magnetic field, A m f p the particle mean free 
path and uu the (ion) collision frequency. This requirement is rarely satisfied in real 
turbulent astrophysical and space plasmas (e.g., in the solar wind, A m f p ~1 AU) and it is 
an observational certainty that turbulence exists at collisionless scales [TT1I2]. The same 
is true in fusion plasmas. Thus, plasma turbulence must be understood in the framework 
of kinetic theory, which evolves the distribution function f s for each species s (= i, e): 

§+^+^(*+^)'§Hi)' < 3 » 

where q s and m s are particle charge and mass, c is the speed of light, the right-hand 
side of ([3]) is the collision integral (quadratic in /), and E and B are the electric and 
magnetic fields, which satisfy Maxwell's equations: 



V ■ E = 4vr ^ Vsris, n s = d 3 v f s , (4) 

u + jcxt) , j = & J d3v v f^ ( 5 ) 



I dE _ An 

c dt c 



f) Pi 

= -cV x E, V • B = 0. (6) 

dt v ; 

In §5§, the external current j ext stands in for the outer-scale energy injection. 

The energy injected into the plasma must be dissipated and converted into particle 
heat. It is in fact a rather subtle issue what this exactly means. Multiplying ([3]) by 
m s v 2 /2 and integrating, we find that the total particle energy satisfies: 

d [d 3 r ^ [ ~ m s v 2 , r d 3 r „ . d f d 3 r E 2 + B 2 

s 

where e = — (1/V) j d 3 r E ■ j cxt is the injected power per unit volume. In deriving the 
above equation, we used Ampere's law Faraday's law (JS)), and integrated by parts 
wherever opportune. Equation ([7D tells us that, unsurprisingly, the change in particle 
energy is equal to the work done on the particles (J d 3 r E ■ j) and that the change in 

where 5u\ is the characteristic relative velocity of fluid elements separated by a distance A and t\ is 
the characteristic nonlinear interaction time [energy- cascade time) at this scale. For a local cascade, 
dimcnsionally, t\ ~ \/5u\ and the Kolmogorov scaling law immediately follows: Su\ ~ (eA) 1 / 3 . 
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the combined energy of the particles and fields is equal to the injected energy. This, 
however, is not yet a statement about heating in the thermodynamic sense of the term 
because the energy exchange described by (I7j) is, in principle, reversible. In order to 
effect irreversible heating, we must change the entropy of the system and that, in a 
closed kinetic system, can only be accomplished by collisions. This result is known as 
Boltzmann's if-theorem [§]: from ([3]), it is readily obtained [36J that the entropy S s of 
species s grows according tc 



dS s d 



dt dt 



d 3 r 



d 3 v f s In f s 



<VV f d*v\nf s (^\ >(). (Si 



V J V dt 



We would now like to assume that the plasma distribution function can be split into 
a slowly changing equilibrium part and a fast changing fluctuating part, f s = F 0s + Sf s , 
that the latter is small, and that its smallness is controlled by some parameter e< 1. 
In the next section, we shall specialize to the case of gyrokinetic turbulence, where 
e ~ uj/Qi, the ratio of the typical fluctuation frequency to the ion cyclotron frequency. 
As we shall see momentarily, the equilibrium quantities can be assumed to vary on a time 
scale ~ (e 2 ^) -1 , much longer than the fluctuation time scale We further assume 
that the collision rate is va ~ ui, i.e., while the dynamics are not collisionally dominated, 
collisions are retained on a par with fluctuations. This can be viewed as a convenient 
ordering prescription on the level of the e expansion [28] and does not prevent one from 
considering the collisional (va ^> uj) and collisionless (ua <C to) regimes as subsidiary 
limits [12]. With these assumptions, (jSJ) implies that the equilibrium distribution is 
a local Maxwellian for each species [HI [36] : F 0s = n 0s (Tiv1 hs )~ 3 ^ 2 exp(— v 2 /v 2 hs ), where 
v t i is = (2T 0s /m s ) 1 / 2 is the thermal speed and T 0s the temperature. For simplicity, we 
shall ignore all spatial gradients of the equilibrium quantities compared to the gradients 
of the fluctuating ones and also assume that the plasma motions are subsonic, i.e., the 
Mach number is small, M = u/v t h s ~ e 1§| 

If we now substitute f s = F 0s + Sf s into (jHJ), use the assumptions explained above, 
and keep only the lowest-order terms in e, we get 



rr, &Ss d 

Tn 



-Os 



dt dt 



V J 2 Js J V J 2F ( 



Os 



d3r f jS-.TosSfs f d5f s \ , f j3 m s v 2 fdF 0s 



V J F 0s \ dt ) c J 2 V dt ) c y ' 

The second term on the right-hand side represents the collisional energy exchange 
between the Maxwellian equilibria of two species and is equal to — n 0s u^' (T 0s — T 0s >), 
where is the appropriate rate of collisions between species s and s' [21]. Equation 
(Q has two key consequences. First, let us average it over times longer than the 
fluctuation time scale but shorter than the equilibrium-variation time scale, uo^ 1 <^ 

4 Boltzmann's function is H = (1/V) J d 3 r J d 3 u /In / = -S, so dH/dt < 0. 

5 These are rarely good assumptions at the outer scale, but, in many astrophysical applications, they 
are increasingly better satisfied as we move deeper into the inertial range [42] . In tokamak plasmas, the 
equilibrium gradients do play an important role, but it is not essential to retain them in the conceptual 
discussion that follows. 
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t <C (e 2 ^) 1 ■ Then the time derivatives of the fluctuating quantities vanish and, noting 
that / d?v (m s v 2 /2)F 0s = (3/2)n 0s T 0s and dn 0s /dt = 0, we get [28] 



I = - / f#v ^ ) - n ,^(T„, - T „), (10) 



dt J V J F 0s V ^ 

where the overline denotes the time average. The first term on the right-hand side is 
positive definite and represents the heating of the equilibrium via collisional dissipation 
of the fluctuating part of the distribution function - - precisely the transfer of the 
fluctuation energy into heat that is the ultimate imperative of turbulence. Note that 
(flOl) is consistent with the ordering assumptions made earlier: the equilibrium evolves 
on the time scale ~ (e 2 co>) _1 , as we have ordered u ss ~ uj. 

The second important consequence of Qj arises if we sum over species and use ([7]) 
to express the first term under the time derivative in ([9]). This gives 



d_ rd 3 r 

dt. V 



^ J 2F 0s 8tt 



The positive definite quantity under the time derivative on the left-hand side, henceforth 
denoted W, will be referred to as generalized energy® Its evolution is determined by 
the competition (or, in a stationary state, balance) of the externally supplied power 
e and collisional dissipation (the negative-definite term on the right-hand side) — the 
latter converts the generalized energy into heat according to (fTUl) . Thus, we have a 
conservation law analogous to (j2j). This suggests a straightforward generalization of 
the view of fluid turbulence outlined in §[TJ to plasma turbulence: its cause and effect 
is the transfer of the generalized energy injected at the outer scale to scales where the 
collisional dissipation can convert it to heat. 

There is, however, an important novel feature here. If the collision frequency is 
small, v ss <C u, the collision term in (1111) can only balance the injected power provided 
the perturbed distribution function develops small-scale structure in velocity space. 
Since the collision operator is a second-order (diffusion) operator in the velocity space, 
we may roughly estimate the smallness of this structure by balancing to ~ u ss v 2 hs d 2 / dv 2 , 
so the correlation scale in velocity space is Sv/v^g ~ (v ss / uj) 1 / 2 . As we shall see in §Hl in 
gyrokinetic turbulence, the emergence of small scales in velocity space is intertwined with 
a cascade to small scales in physical space. Thus, in the same way as fluid turbulence 
could be described as the energy cascade, plasma turbulence is a cascade of generalized 

6 We use this term to emphasize the role of W as the cascaded quantity in plasma turbulence (see 
below). The importance of its conservation for plasma turbulence was realized by several authors 
[TBI [22, 52, 28, 43J, who refer to it as the "generalized grand canonical potential" or free energy. The 
latter term is perhaps physically the most appropriate because it flags the interpretation of W as the 
work content of the particles + fields system. The part of W that involves 5f s is equal to — Toa5S s , 
where 6S S is the perturbed entropy. In an exactly collisionless plasma, ((9]) and ([7|) show that any work 
done on the plasma simply increases this quantity. Any increase in the "equilibrium" entropy that 
might appear to be heating is then, in fact, compensated by a decrease in the perturbed entropy, so 
this "heating" is reversible. See [2U [33l US] for discussion of the entropy production in plasmas. 
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energy, or a kinetic cascade — this cascade occurs in phase space, reaching towards 
small scales both in physical space and in velocity space. 

If the heating is always ultimately collisional, what then is the status of the 
collisionless (Landau) damping [35] as a dissipation mechanism for (homogeneous) 
plasma turbulence? Collisionless damping does not appear explicitly in (TTT]) because 
what it does is, in fact, redistribute the generalized energy: the energy of electromagnetic 
fluctuations (E 2 + B 2 ) is converted into entropy fluctuations (T Qs 5f 2 /2F 0s ). In order 
for any actual heating to occur (i.e., for the fluctuation energy to be lost irreversibly), 
this perturbed entropy has to be transferred through phase space to collisional scales. 
There are two ways in which this can be accomplished: linear and nonlinear. The first is 
the well known [231 IE] phase-mixing mechanism associated with the so-called ballistic 
response in the perturbed distribution function: the linearized kinetic equation ([3]) has 
the homogeneous solution Sf s oc e~ lk ' vt [35], for which d8f s /dv ~ kt5f s , i.e., there is 
a secular growth of the velocity-space derivatives and the collisions become important 
after a time t ~ {kvth s )^ l {uj / 'u ss ) 1 ^ 2 . In fact, as anticipated in [U] and as we will show 
in §Hl the linear phase mixing can be superceded by a faster nonlinear mechanism that 
cascades the perturbed entropy to collisional velocity scales over times t ~ 

Finally, we note that one can make a plausible argument in favour of an effectively 
irreversible "collisionless heating" in the sense that the distribution function may become 
so convoluted in phase space that it is effectively impossible to unscramble it and 
the entropy of an approppriately defined "coarse-grained" distribution is increased. 
Discussions of this process and the difference between such effective irreversibility and 
the exact irreversibility for which collisions are necessary have continued since the 



salient physical fact that until the collsions can act, the negative entropy necessary to 
compensate for the increase in the coarse-grained entropy is stored in the fluctuations 
of the perturbed distribution function and that these fluctuations are explicitly present 
in the overall generalized energy budget (fTTj) — a conservation law that underpins the 
interpretation of plasma turbulence proposed here. 

3. Gyrokinetics and the Many Forms of the Kinetic Cascade 

Our treatment so far has not been specific to gyrokinetic turbulence. However, the 
particular mechanism of kinetic cascade in phase space we intend to discuss in §H] will 
be. Thus, we now briefly introduce the gyrokinetic approximation and describe the 
forms the kinetic cascade from macro to microscales takes in gyrokinetic turbulence. 

It is nature's gift to plasma physicists that magnetized plasma turbulence both in 
fusion devices and in space appears to consist mostly of fluctuations whose frequencies 
are much lower than the ion cyclotron frequency, uj -C f2j, even as their spatial scales 
perpendicular to the magnetic field can be as small as or smaller than the ion gyroscale 

7 Understanding the role of phase mixing and weak collisions in converting wave energy into heat is 
practically important, e.g., in the theory of RF heating; see, e.g., [47l [6] and references therein. 



birth of quasilinear theory to the present 
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Pi = i>thi/^V This low-frequency character of the turbulence is intimately related to the 
tendency of plasma fluctuations in a dynamically strong magnetic field to be spatially 
anisotropic, with fcii <C k±. Let us briefly explain why. 

The structure of plasma turbulence is set by the interplay of parallel 
linear propagation effects (waves, particle streaming) and perpendicular nonlinear 
decorrelation (turbulent cascade). It is crucial to understand that, while anisotropic, this 
is an essentially three-dimensional situation. For fluctuations with a given perpendicular 
correlation length, the parallel correlation length is set by the distance a wave (or 
streaming particles) can travel during one perpendicular correlation timejf] A good 
example of this principle is the Alfvenic MHD turbulence, where it is known as the 
critical balance [20l |2T] . Alfvenic turbulence is the predominant type of turbulence in 
finite-beta plasmas at scales above the ion gyroscale (the "inertial range") irrespective 
of the degree of collisionality — this statement can be proven analytically [SJ H2] and 
there is ample evidence in its favour from measurements in the solar wind [TTJ [2]. 
Alfvenic fluctuations have velocities and perturbed magnetic fields u± ~ 5B±/y/Anminoi 
perpendicular to the mean field B = B z. Their decorrelation rate is ~ k±u±, while the 
characteristic propagation frequency is u = k\\VA, where va = Bq/ ^f^Km^n^. In critical 
balance, k\\VA ~ k±u±, so ku/k± ~ u±/va <C 1. If the Alfvenic cascade from the outer 
scale to the ion gyroscale respects this principl^] (and there is numerical [12j [37] and 
observational [25] evidence that it does), the fluctuation frequency at k±pi ~ 1 will still 
be low compared to the ion cyclotron frequency: u/Qi ~ fciif^/f^ ~ (k\\/k±)k±pi/ \/A ^ 
1 (we assume moderate values of /?;). 

The gyrokinetic approximation can now be constructed by using the critical balance 
explicitly as the ordering prescription: e ~ k\\/k± ~ uj/VLi ~ u±/va ~ q s ip/T 0s ~ 
5B±/B ~ 5Bu/B ~ 8f s /F 0s , where <p is the scalar potential^ The Vlasov- 
Maxwell equations ([3])-([6]) are expanded in e and averaged over the particle gyromotion 
[T7| 1281 [TO] . As a result of this procedure, the perturbed distribution function splits 
into the Boltzmann response and the perturbed distribution of particle gyrocentres: 
Sf s = —q s (pF 0s /T 0s + h s (t, R s , v ±, v\\), where R s = r + v± x z/Q s is the gyrocentre 

8 Clearly, perpendicular planes separated by longer distances cannot remain correlated, which rules out 
the two-dimensional limit (linear frequency <C nonlinear decorrelation rate). Decorrelation at shorter 
distances gives rise to weak turbulence (linear frequency 3> nonlinear decorrelation rate), which tends 
to produce a cascade towards smaller perpendicular scales, where the nonlinear decorrelation rate again 
becomes comparable to the wave frequency [HJ IIH] • 

9 A one-paragraph review of the Goldreich-Sridhar 1995 MHD turbulence theory [201 121) : Making 
the same assumptions as in Kolmogorov's theory (footnote [3]) except isotropy, we have, for Alfvenic 
velocities, 8u\/t\ ~ e, where A is now the perpendicular scale. If the critical balance holds, 
v a/1\\\ ~ 8u\/\, where is the parallel correlation length of these fluctuations. Since this means 
that only one time scale is present in the problem, we must have t\ ~ X/Su\ and thus recover the 
Kolmogorov scaling: 5u\ ~ (eA) 1 / 3 . Using this and the critical balance, we find the scaling relationship 
between the perpendicular and parallel scales: lt\x ~ /q^ 3 A 2 / 3 , where Iq = v A /e (but see [7] for a 
version of this theory giving rise to different scalings). Thus, there is a cascade both in the parallel and 
perpendicular directions, but the aspect ratio increases as we move deeper into the inertial range. 
10 The Alfvenic velocity perturbation is the E x Bq flow: u± = z x V ±cip/ Bq. 
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position. In a uniform magnetic field B , the gyrokinetic equation for h s is 

dh s . dh s c f/ . q s F 0s d(x)R s ( dh s \ 

+ v \\ -zr + 77" {(x)R s ,h s } = — 77—+ Hrr , (12) 



where \ = V ~ v • A/c, B = B z + 5B, 5B = V x A, V • A = 0, and (■ ■ -) Rs 
is the gyroangle average at constant R s . The vector potential A is recovered from 
([5J) neglecting the displacement current. The scalar potential y is found from the 
quasineutrality condition: neglecting V • E in (j3j) and separating the Boltzmann 
response, we have 

Ef^ n o s =y> / dMM., as) 

where (• • -) r denotes gyroaveraging at constant r (the velocity integral is at constant r). 

Gyrokinetics helps make the problem of kinetic cascade numerically [301 EQ] arid, in 
certain limits, analytically [121139] tractable because all high-frequency physics (u > fl*) 
is systematically ordered out and the gyroaveraging reduces the phase space from 6D to 
5D. However, it is still a fully kinetic system and everything that was said about heating 
and the kinetic cascade in §[2] remains valid. The generalized energy conservation law 
( ITTT) for gyrokinetics takes the following form [28l [42] : 



dW d f d 3 r 



dt dt 



r&r y( f d s v TosQ§r 9 5 V^ A , \SB\ 
J V \J 2F 0s 2T 0s J 

f ,3 f d 3 R s To s h s ( dh R 



in 



As we explained in §[21 the generalized energy injected at the outer scale has to be 
transferred (cascaded) through phase space eventually to reach the collisional scales. If 
we conjecture that this transfer is local in scale space, we can ask what forms the kinetic 
cascade takes in several distinct physical regimes separated by the characteristic plasma 
microscales: the mean free path, the ion and the electron gyroscales. It turns out that 
in each of the asymptotic limits k\\X m fp 1, k\\\ m f P 1, k±pi <C 1, k±pi ^> 1, etc., 
the kinetic cascade splits into several non-energy-exchanging channels corresponding to 
cascades of distinct plasma fluctuation modes, some of which are familiar from fluid 
models of plasma turbulence and some are new. As the characteristic scales are crossed 
(k\\\ m f p ~ 1, k±pi ~ 1), these channels join together into a single cascade, which then 
splits again, but in a different way, as another asymptotic limit is reached. All these 
asymptotic limits are worked out in detail in [12]. Here we briefly summarize their 
role as a route for the generalized energy to reach the ion gyroscale (at which point 
interesting things start happening in the phase space). 

Let us imagine that energy is injected at scales larger than both the mean free path 
and the ion gyroradius. As the cascade takes the energy to smaller scales, anisotropy 
and critical balance are established, so the gyrokinetic approximation applies [29l [27]. 
In the inertial range (k±pi <C 1), the energy cascade is split into two main channels: the 
Alfvenic turbulence (5B±, u±), which is described by the Reduced MHD equations [IE] 
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regardless of the collisionality [4"Tl |4"2"] and the "compressive" component {5n, un, SBn). 
The Alfvenic cascade is split into two cascades corresponding to the two directions of 
propagation of the Alfven waves. While the nonlinear interaction is always between 
the "+" and "— " waves, it is of "scatter" type, so no energy is exchanged between the 
two cascades. The compressive fluctuations are passively mixed by the Alfven waves, 
again without energy exchange. In the collisional limit (fc||A m f p <C 1), the compressive 
cascade is split into three channels: the "+" and "— " slow waves and the entropy-mode. 
As feiiAmfp ~ 1 is approached, these three are mixed together and remain mixed for 
k\\\ m fp 3> 1. Dissipation and collisional heating can occur at this transition because in 
the fluid limit, the collisional term in (1111) can be activated by small deviations of the 
distribution function from a Maxwellian — the smallness of the collision rate in this case 
is overcome not by a velocity-space cascade but by the fact that the non-Maxwellian 
part of the perturbed distribution function is proportional to k\\/vu. At collisionless 
scales, the compressive fluctuations experience the Barnes (transit-time) [3j H7] version 
of Landau damping — as discussed in §EJ this transfers the energy associated with the 
compressive fluctuations into ion entropy fluctuations. 

As the inertial-range cascade transfers energy to scales around k±pi ~ 1, its Alfvenic 
and compressive components cease to be decoupled from each other and all fluctuations 
are subject to Landau damping (see [28] for details of linear gyrokinetics). What emerges 
on the other side of this transition, at k±pi ^> 10 is a cascade of generalized energy again 
split into two channels: the fluctuations polarized as kinetic Alfven waves (KAW) (they 
satisfy fluid-like equations closely related to Electron MHD [3UI12]) an d, energetically 
decoupled from them, the ion entropy fluctuations (T 0i hf /2F 0i ). The latter carry the 
part of the inertial-range energy that was Landau-damped at the ion gyroscale and, 
possibly, also in the inertial range (for the compressive fluctuations). How it becomes 
ion heating is the subject of §HJ The KAW cascade takes the energy to the electron 
gyroscale, k±p e ~ 1, where it is converted by Landau damping into electron entropy 
fluctuations (T Qe hl/2F 0e ), eventually giving rise to electron heating in a way analogous 
to the ion case discussed below. 

4. Nonlinear Perpendicular Phase Mixing and the Entropy Cascade 

In order to introduce the concept of the phase-space cascade of entropy in the simplest 
possible setting, we will consider the extreme case where all of the fluctuation energy 
arriving to the ion gyroscale from the inertial range is converted into entropy fluctuations 
by Landau damping, i.e., we will neglect the KAW component of the dissipation-range 
turbulence £f] Furthermore, we will use the Boltzmann-electrons approximation, which 

11 In space physics, this is called the "dissipation range," a historical misnomer dating back to the times 
when it was not appreciated that it can contain dissipationless cascades. 

12 In the presence of KAW turbulence, the entropy fluctuations are passively mixed by KAW. This case 
can be treated in a way analogous to what we do here [42] . At moderate values oi Pi, a KAW cascade 
is probably a good description of the dissipation-range turbulence in the solar wind [29J, [30] . The case 
without KAW may be more relevant at low and high fii because ion Landau damping is quite strong 
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Figure 1. The nonlinear perpendicular phase mixing: the gyrocentre distribution 
function at a given point Ri is mixed in a decorrelated way by E x B flows 
gyroaveraged over ion orbits whose radii an d Uj_/f2j) differ by more than the 

flows' correlation length. 



is justified to the lowest order in the mass-ratio expansion as long as we stay above 
the electron gyroscale, k±p e <C 1 [151 02J. These approximations mean that we have 
5f e = eipF 0e /T 0e , i.e., h e = 0, while hi satisfies the electrostatic version of ([12]) (x = V 9 )- 
The resulting system of equations follows from (fl2l) and (TT3|) : 



where Z = qt/e, r = T^jT^ and hi(k) is the Fourier transform of hi(Ri). 

As we explained in §[2] in order for the collision term in ([15]) to become non- 
negligible, small-scale structure has to be generated in the velocity space with Sv/v t u ~ 
(z/jj/u;) 1 / 2 . One route to such small scales is via the parallel (linear) phase mixing, whose 
role in plasma turbulence has been well established for some time [23] EH E21 E2]: the 
ballistic response hi cx e Hi* gives rise to secularly growing gradients dhi/dvu ~ kvthi 
and, therefore, small scales in parallel velocities: Sv\\ ~ l/k\\t. 

The other, perpendicular, phase mixing mechanism is nonlinear [TJJ02]. 111 (TT5|) . 
the nonlinear term represents mixing of the ion distribution in the gyrocentre space by 
the gyroaveraged ExB flows. Like any random mixing, this produces small scales in Ri. 
It also produces small scales in Vj_ for the following reason. Consider ffT5l) taken at two 
different values of velocity, v± and v' ± . Then hi(Ri,v±) and hi(Ri,v'j_) will be spatially 
mixed by the gyroaveraged ExB velocity field given by (<^)K i (f±) and (</?) ^(f^) , 
respectively. These gyroaverages come from spatially decorrelated fluctuations of ip if 
the difference between the gyroradii v±/Qi and v' ± /Qi is larger than the perpendicular 
correlation length l/k± of <p (figured]). If this condition is satisfied, hi(Rt,v±) and 
hi(Ri,v' ± ) are mixed by decorrelated fields and are, therefore, themselves decorrelated. 

then 28 . It is also interesting in the context of electrostatic microturbulence (ITG. ETG, drift waves) 
prevalent in fusion plasmas [TJl [T5l [26] . 




(15) 



V ZJ T 0i n 0i J n 0i J \ Ik J 



(16) 
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Thus, small-scale structure of ip in the physical space gives rise to small-scale structure 
of hi in the velocity spaced the correlation scale in the velocity space is 



Sv i 1 



Vthi Pi 



' < 1 when k ± pi > 1. (17) 



k ± p 

The small-scale structure of hi in the gyrocenter space gives rise to similarly 
small-scale structure of ip in the physical space. Using ( fl6l) . they can be related as 
follows. For k±pi 3> 1, the Bessel function in the velocity integral is Jo{k±v±/Qi) ~ 
(2{l i /nk±v±) 1 / 2 cos(kj_v_L/^li — tt/4), i.e., it oscillates in v±_ with the period Sv±/v tili = 
27r/k±pi. But, according to (j!7p . this is also the correlation scale of hi(k) in velocity 
space. Assuming that the velocity integral accumulates as a random walk and taking 
into account also the factor of 1/ \Jk^pi coming from the Bessel function, we have 

Zeip(k) v 3 hi 1 f$v x \ 1/2 ,,,. v^hi(k) 



, , , hAk) 

Toi n 0i y/k±pi V «thi / n 0i k ± pi 

The gyroaveraged potential is then Ze{ip) ^(k) /T 0i = ZeJ (k±v±/Qi)ip(k)/T 0i ~ 
( v thi/ n oi)hi(k) I (k±pi) 3 / 2 , and so the perpendicular mixing of the particle distribution is 
a fully nonlinear process. 

This process can be understood as a local (in scale) ion entropy cascade and a 
Kolmogorov-style scaling theory can be constructed for it. Recall that the gyrokinetic 
equation (1151) has a conservation law given by (114j) with h e = and SB = 0. In 
view of f|T8|) . Z 2 e 2 ip 2 n 0i /2TQi J d 3 v T Qi h 2 /2F Qi) so the entropy of the perturbed ion 
distribution is conserved individually. This is, in fact, obvious also from ( fl5l) : again 
using ([TBI) , the inhomogeneous term on the right-hand side is negligible for kj_pi ^> 1 
and J d 3 Ri h 2 is clearly a conserved quantity but for collisions. Denoting by ip\ and hi\ 
the characteristic fluctuation amplitudes at some perpendicular scale A <C pi and by t x 
the corresponding cascade time, we may write (cf. footnotes [3] and [9]) 

i^&~e, TA ^y /2 A2 „ ft /2 * i/2 "oi (19) 

n i r x \\J cipx/B v* hi h iX 



where we used ( fT81 to get ap\j B ~ v^ hi hix\/n 0i . Combining these relations, we finco 

Ze^^J^ TX ^'IA^!1 (20 ) 

iA „3 ,1/3 ' rp 5/6,1/3' A „ , . ' ^ 

where Zo = m i n oiVthi/ £ - These scalings correspond to a /c_ l _ 4/ ' 3 spectrum of hi and a fc_|_ 10 ^ 3 
spectrum of ip. Encouragingly, these predictions seem to be corroborated by numerical 
simulations of electrostatic gyrokinetic turbulence in two spatial dimensions |50j . 

Now let us revisit the question of parallel phase mixing. In our discussion of the 
perpendicular cascade so far, we have ignored the presence of the parallel propagation 
(particle streaming) term in ( ITBI) . In a formally 2D situation, i.e., when uj ~ r^ 1 ^> k\\V\\, 

13 Note that this nonlinear perpendicular phase mixing mechanism was first recognized in [14] : in their 
gyrofluid closure formalism, it manifested itself as the growth of high-order v±_ moments of hi . 
14 It is also possible to derive exact scaling results analogous to Kolmogorov's 4/5 law, which prove to 
be consistent with (J20J) [39]. 
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this is, of course, allowed and the perpendicular scalings derived above should hold. 
However, it is more likely that the parallel scale of the fluctuations will adjust to 
their perpendicular scale according to the critical balance principle explained in §[3j 
for the fluctuations with perpendicular scale A will be such that particles can stream 
across the distance l\\\ in one nonlinear decorrelation time T\. Using ( 1201) . this implies 
h\ ~ ^thi r A ~ ^o^P^ 3 ^ 1 ^ 3 (cf- footnote [9]). As we explained at the beginning of 
this section, the typical parallel correlation scale in velocity space produced by the 
parallel phase mixing is 8v\\ ~ l/k\\t. Therefore, after one perpendicular cascade 
time, no appreciable refinement of the parallel velocity-space structure is achieved: 
8v\\/vtbi ~ l\\\/vthiT\ ~ 1. In contrast, in v±, one cascade time is enough for the entire 
cascade down to the collisional cutoff (to be calculated in §|5]) to be set up. Thus, the 
linear parallel phase mixing is much less efficient than the nonlinear perpendicular one. 



5. Conclusion: Dissipation Achieved 

Let us now come back to the original motivation for the above developments: the 
necessity to understand how the distribution function is brought to collisional scales in 
the velocity space. We have seen that this is done by transferring the energy injected 
at the outer scale down to the ion and electron gyroscales via a multichannel cascade 
of generalized energy through phase space. Below the ion gyroscale, the phase-space 
nature of the cascade becomes particularly manifest as the ion distribution function 
simultaneously develops small scales in the gyrocentre and velocity space via a nonlinear 
perpendicular phase mixing process. We have described this process as a Kolmogorov- 
like turbulent cascade enabled by a constant flux of ion entropy and derived scaling 
relations for the fluctuations of the distribution function and the electric potential. 

Using these scalings ( 1201) . let us now estimate the collisional cutoff in phase space. 
As we explained in §|2l the collisional scale is reached if the velocity-space correlation 
scale is Sv/v t u ~ (vu/lu) 1 / 2 . Using (1T7I) and estimating u ~ , we get 

x 1 7 1 / 5 2 / 5 

^ ~ -J- ~ (vt ) 3/5 ~ l ° P * (21) 

where Iq = minoiV 3 hi /e and r Pi ~ (min 0i p 2 / e) 1 ^ 3 is the fluctuation time scale at k± P i ~ 1. 
This formula is perhaps our most consequential result for numerical applications: 
it tells us what it means to have a well-resolved gyrokinetic simulation of plasma 
turbulence and shows that the resolution requirements in the gyrocentre and velocity 
spaces are fundamentally linked. In this context, it is clear that adequate modeling of 
collisions [TJ [5] and controlled velocity-space resolution [1] are imperative for gyrokinetic 
simulations. 

No matter how small the collisional cutoff (12T1) is, all of the energy channelled 
into the sub-gyroscale entropy cascade will reach this cutoff in finite time — roughly 
the nonlinear interaction scale r Pi evaluated at the ion gyroscale. Since the process is 
nonlinear, this time is amplitude dependent. If the principle of critical balance (§[3]) 
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holds at the ion gyroscale, t Pi should be roughly equal to the linear parallel propagation 
time scale at k±pi ~ 1. Importantly, the time to reach the collisional cutoff does not 
depend on the collision rate — just like in hydrodynamic turbulence (§[[]), the time to 
reach the viscous scale is the turnover time at the outer scale, independent of viscosity. 

Another important conclusion is that the dissipation range (k±pi > 1), even in the 
absence of kinetic Alfven waves, is filled with electrostatic fluctuations due to the ion 
entropy cascade. This is a purely kinetic effect invisible in any fluid models. In fusion 
plasmas, this may be relevant for identifying the nature of electrostatic fluctuations 
found between the ion and electron gyroscalesf^l In space physics, the great variability 
of the observed spectra in the dissipation range [44] might be speculatively attributed 
to varying proportions of energy contained in the entropy and kinetic-Alfven-wave 
cascades [42]. 

These results are only the first glimpse of what one finds if one adopts the view of 
plasma turbulence as a kinetic cascade in phase space. We believe that further studies 
conducted in this vein, both numerical [30], [50] and analytical [l2l [39] , will unveil much 
new physics and many new and tantalizing questions. 
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